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An essentially analytical theory of atmospheric cosmic -ray prop- 
agation is developed on the basis of a phenomenological model of 
hadron-nucleus collisions. This model correctly predicts the sea- 
level cosmic-ray nucleon, pion and muon spectra, the cosmic-ray 
ionization profile in the atmosphere, and neutron flux and density 
profiles in the atmosphere. It is concluded that the large scale 
properties of atmospheric cosmic-rays can be accurately predicted on 
the basis of a purely nucleonic cascade as a result of which all 
secondaries are mediated by pion production. 

Implications for energy independence of cross sections, the 
recent 70 GeV results from Serpukhov, and nucleonic relaxation rates 


in the atmosphere are discussed. 


INTRODUCTION 

This paper attempts to establish the 
physics on which the large-scale, time- 
averaged, one -dimensional properties of 
galactic cosmic-rays in the atmosphere 
depend. 

The point of departure for this 
theory is a phenomenological model of 
high energy nucleon-nucleus collisions 
which can be applied to analytical trans- 
port theory. Particle spectra, fluxes, 
densities and ionization calculated from 
the theory yield good agreement with 
measured values indicating the adequacy of 
the nuclear model and of the supporting 
cosmic -ray and geophysical data. All 
compar isons are on an absolute basis. 

Preliminary results of this work have 
already been reported comparing calcula- 
tions and measurements of various 
components of cosmic-ray ionization in the 
lower atmosphere (<2.5 km elevation) at a 
geomagnetic latitude of 51° (ref. 1) . 


ATMOSPHERIC MODEL 

The atmosphere is assumed to be a 
flat slab 103 3 g/cm* 3 thick with a constant 
scale height of 6.7 km. It is assumed to 
be composed of a single nuclear species 
with an atomic weight of 14.48, atomic 
number of 7.31, and an ionization poten- 
tial of 86.8 volts. Because oxygen and 
nitrogen are so close in the periodic 
table, this simple assumption yields the 
correct nuclear data. The density of the 
atmosphere is 

P = r/H (1) 

where 

p is the dehsity in g/cm 3 , 

r is the depth in the atmosphere in 
g/cm 8 , and 

H is the scale Height in cm. 

Since the mean free path for decay of 
a charged particle is 
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is the mean free path for decay in 
g/cm 3 of a particle of type q, 
is the mass of the particle in 
MeV/c s , 

is the momentum in MeV/c, 
is the velocity of light, and 
is the mean life in the rest frame , 
in seconds. 


we have the useful result that 

= ^q cr /Cqj (3) 

where 

C q = mqC S H/cT g . 

THEORY OF THE ATMOSPHERIC NUCLEONIC 
CASCADE 


The Boltzman Equations 

In this paper, atmospheric cosmic- 
ray fluxes will be obtained as analytic 
solutions to an approximate form of the 
Boltzman equations. The Boltzman equa- 
tions for the nucleonic cascade in the 
atmosphere are 

B q cpq(r,E,Q) = S g j , (4.1) 

(q = P>n,nsj = p,n), 

B^cp^(r,E,n) = S^i, 

Va (r ’ E ’^ = S af3> 

(a = y,e:f3 = y,e,n ,n) , 


Bg = Q 


• V + CT a + — 

g V r 

(q = P,n,TT,n,eY) 



(4.2) 


S qj J 4n d ^ J E dE B ff qj p qj ( e B - E,n'*0 # ). 

•cpj (r,E B ,n') (4.3) 

(qj = np, pn, nn, rrp, pin* ept, yn 0 , ey, ye) , 

a n = C p = c n = C e = C Y = *n = k n° = k y = °' 
where 

r is the depth in the atmosphere in 
g/cm s , 

E is the particle kinetic energy in 
_ MeV, 

0 is the unit vector in the direction 
of particle travel, 

cpg is the particle flux of type q per 
second per steradian at a depth 
r with a direction 0, 


ct is the total cross section for 
* absorption of a particle of type 
q in cm a /g, 

is the cross section for the pro- 
duction of particles of type q 
from collisions with, or decay 
by, particles of type j in 
cm*/g, 

k is the stopping power of a charged 
particle of type q in air, in 
MeV cm a /g, and 

F . is the number of particles per MeV 
3 per second per steradian at E 

and Cl resulting from a collision 
with or decay by a particle of 
type j at E fi and Q 7 . 

The subscript tt implies application 
to all pions, charged and neutral. The 
subscript tt* implies application to the 
charged pions, and rr° to neutral pions 
only. 

Nucleon-Nucleus Collisions 

In eq. (4) it has been assumed that 
atmospheric cosmic -rays propagate by means 
of the nucleonic cascade in an exponential 
atmosphere [see Fig. 1 of O'Brien (ref. 1)]. 
Thus all other secondaries result from 
nucleon-nucleus collisions. The following 
reactions are considered. 


p + air - VpP + v n n + v n ± tt* + v n o tt 
n + air -* VpP + v n n + v n ± tt* + v n ° tt 
tt* -» \i + v 

tt° -► 2y -> electromagnetic showers 
p. -* e + 2v -* electromagnetic showers 


o 

o 


(5) 


where Vj are the multiplicities of j type 
resulting from the collision of a nucleon 
with a nucleus of air. The influence of 
kaon production has earlier been shown to 
have a small influence on ionization, and 
it is neglected here, at a considerable 
saving of computer time (ref. 1) . 


The nucleon-nucleus reactions of eq. 
(5) will be considered at high energies 
only, because the mathematical form of the 
approximation to the Boltzman equation to 
be obtained is only applicable at Energies 
above about 0.1 GeV*. 
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It is required (a foreseen mathemat- 
ical result motivates the choice of 
function) that 

e b 

s qJ ■ ? D(E » ' v 

r rr 

Ggj = 2tt j o d9 sine F q j (E B -» E, 
o = cose 


where 

Ggj is the secondary production 

spectrum of type q particles 
integrated over the solid angle, 
Iq is an arbitrary constant depend- 
ing on q, 

i and n are arbitrary constants 

which must be the same for all 
j,q, and 

U(x) is the Heaviside function 

[u (x < 0) =0, U (x £ 0) = 1], 
r\ is a lower energy limit below 
q which secondary particle 
production is cutoff. 

The formula for G q j is certainly very 
crude, however it is suitable to represent 
the behavior of the partial inelasticities 
and multiplicities associated with high 
energy nucleon-nucleus collisions. 
Inelasticities are known to vary quite 
slowly with energy, and to become essenti- 
ally constant at energies of a few 10's of 
GeV (ref. 2) . If n = it 1, then G q j can 
be rewritten in terms of a constant 
partial inelasticity, K_, for the produc- 
tion of a type q particle 

E X 

Ggj = (1 - -e)Kq ^ U(E B - T^) 

Hagedorn and Ranft (ref. 3) have 
calculated Kg using the statistical model 
for p-p collisions at 12.5, 18.8, 30 and 
300 GeV/c and it would be quite convenient 
to use these results for Kg. However, 
Alsmiller and Barish (ref. 4) have shown 
that the secondary production spectrum. 


*Low energy nucleon transport (<0.1 GeV) 
is chiefly neutron transport. Low energy 
neutron and electromagnetic shower trans- 
port are treated only very roughly here 
due to the limitations of the transport 
theory to be described. It is intended in 
the near future to apply the S n method to 
this problem, and treat it much more gen- 
erally. The analytic theory presented 
here is if less general, quite accurate at 
high energies, quite simple and rather 
transparent. 


Gqj, softens with increasing atomic 
weight as a result of the intranuclear 
cascade. This effect has been simulated 
in eq. (7) by making Kg a function of 
atomic weight (ref. 5) . This is, unfortu- 
nately, at the expense of the conservation 
of energy in high atomic weight nuclei, 
although for air, energy is conserved 
reasonably well. The intranuclear cascade 
conserves energy as G qj softens, by the 
emission of low energy particles but these 

low energy particles are not important to 
deep penetration calculations at high 
energies. Values of Kg interpolated from 
among the values obtained earlier are 
given in Table I (ref. 5) . 


Table I. - Partial inelasticities for 
proton-air collisions. 


q 


p 

0. 211 

n 

0. 211 

TT + 

0. 180 

TT" 

0. 112 

TT° 

0. 180 

K + 

0. 034 

K" 

0.022 

K° 

0. 034 


Experimental values of partial in- 
elasticities in the energy range 0.1 to 20 
TeV have been obtained in nuclear emulsion 
from the Brawley and I.C.E.F. emulsion 
stacks (ref. 2) . The quantities measured 
were? K^, the energy that goes into new 
charged particles, and hence 

Kch = K n + + IV + K k + + %-? (8) 

Ky is the energy that goes into photon 
production, and so 

*y = KjtO; (9) 

Kb is the energy that goes into long lived 
neutral particles. It is assumed that 
these are neutral kaons. This is a small 
number, and any error that may arise here 
is un inport ant to the study. The energies 
involved in the measurement are high com- 
pared to, the largest rest masses involved 
and thus they may be neglected. 


the 

( 6 ) 

0) 
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In Table II, the predictions of 
O'Brien (ref. 5) (which in conjunction 
with eq. (7) will be referred to as the 
power law model hereafter) is compared 
with these data, and with the prediction^ 
of other nuclear models. 

These other nuclear models have all 
been compared against, and in some cases 
based upon, accelerator target yields, 
mostly at small angles and at high 
secondary momenta. The well known CKP 
model (ref. 6), the Trilling model (ref. 

7) , the extrapolation model (ref. 8) , and 
the TRB model (ref. 9) are considered. 

Ran ft and Borak (ref. 9) have modified the 
formulae used by Trilling (ref. 7), and 
this is referred to as the TRB model. 


Table II. - Partial inelaaticiti.es at very high energies. 



Exper intent al^ 
(0. 1-20 TeV) 

Statistical® 
model 
(300 GeV) 

Power Loir/ 7 
model 

M 

TRET 
model 
(20 GeV) 

Bxt r ap ol at i on - ® 
model 
(200 GeV) 

Trilling® 

model 

i=J 

CKP® 

model 

K ch 

0. 31* . Ob 

0. 

0.303 

0. 23 

0. Ifab 

U. 21 

0.36 

*V 

0.16 

0.186 

0.157 

0.14 

0.084 

0. 09 

- 

Ko 

0. 03 

0.035 

0.029 

* 

* 

0.01 

- 

Total 

0. 50±0.07 

0.585 

0.469 

- 

- 

0. 31 



a nuclear emulsion. 

B hydrogen. 

a air, kaons neglected. 

£ aluminum, kaona neglected. 


The considerable variation is occa- 
sion for surprise. The statistical model, 
the power law model which is an adaption of 
it, and the CKP model are in agreement with 
the cosmic-ray emulsion measurements, 
though all the models appear to agree 
within a factor of two. This is probably 
a consequence of the fact that most of the 
accelerator target data, and much of the 
physical interest, is at small forward 
angles and at large secondary momenta, and 
this does not determine Kq precisely enou^i. 

The secondary particle multiplicities 
of the power law model depend on the value 
chosen for i . Since eq. (7) has an "infra- 
red" divergence, it it not suitable for 
total particle yields. However, shower 
particles in an emulsion produced by high 
energy nucleons have finite lower energy 
limits of 80 MeV for mesons and 500 MeV 


for protons (ref. 10) . The power law model 
yields as the multiplicity of particles 
above a lower limit F 

Vq(E B ) = ^(1 - /je]Kg[ (Eg/F) ^ - 1]. 

Using the data of Meyer et al. (ref. 
11) , the best value of l, in the least 
squares sense, was chosen. The procedure 
is described in a somewhat more expanded 
way in O'Brien (ref. 5). In Table III, 
the predictions of shower particle produc- 
tion by various nuclear models are given. 
The power law result is not really a 
prediction but a fit. Since the statisti- 
cal model and the extrapolation model 
cannot be manipulated without the proper 
computer codes, their reported total 
charged particle production is given in 
place of the shower particle multiplicities. 
As the Trilling model does not give back 
emitted particles correctly (ref. 7) it 
has been omitted from consideration here. 

Some of the differences in Table III 
are certainly due to differences in the 
target nucleus and in the lower energy 
limit. But it is clear that the power law 
model agrees with the experimental data as 
well as the other models. 

Multiplicities and inelasticities are 
averaged quantities related to hadron- 
nucleus collisions. Matters are different 
when one considers the form of Gqj predic- 
ted by the various models. Figure 1 
exhibits G n ±p(E B , E) for Eg equal to 10 
and 100 GeV protons incident on air 
calculated using the power law and TRB 
prescriptions, and the CKP prescription 
for protons on hydrogen. The CKP and TRB 
models agree at high secondary momenta, 
but differ elsewhere. The greater 
sophistication of the TRB model can be 
seen in the graph. For instance, the 
inflection point at G^ipflOO, 18) corre- 
sponds to the transition from energetic 
pions resulting from isobar decay to low 
energy pions emitted isotropically in the 
center of mass. The power law model 
overestimates pion production at high 
secondary momenta but underestimates at 
low secondary momenta. 

In Fig. 2, secondary proton produc- 
tion predictions Gp p (E B , E) are exhibited 
for protons on air, again for Eg equal to 
10 and 100 GeV. The relative crudity of 
the power law model clear. It under- 
estimates proton production at high 
secondary momenta and overestimates at 
low secondary momenta. It is probably 
more significant for the cascade 
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APPROXIMATE SOLUTION TO THE 
BOLTZMAN EQUATIONS 



Fig. 1. Three models of charged pion production 
spectra from proton-air collisions. 


Hadrons and Muons 

A solution will be obtained for a 
sort of "Green's function", that is for 
incident nucleons homogeneous in energy 
and angle, of unit strength per steradian, 
the integral of which over the cosmic-ray 
primaries yields atmospheric cosmic-ray 
fluxes per steradian per second. 


Making the "straight-ahead" 
approx imat ion 

(1 - O' • 0) 


F qj " G qj 6 


2rr 


and eq. (4.3) becomes for the hadron 
component 


- j!° 


dEg c Gqj (Eg, E)cpj (r. Eg, 


(qj = np, pn, nn, np) 


(ID 


0) 

( 12 ) 



Figure 2. Two models of proton production 
spectra from proton-air coll ij ions. 


where cr is the geometric reaction cross 
section in cm 2 /g, and for the muons 
(ref. 13) 

S |xn ± = 2n~^ p 7 r cp TT ± ^ r(m rT ±//l V ) e B’ 

(13) 

where 0 is the cosine of the zenith angle 
of the incident radiation and Eo is the 
energy of the incident nucleon. 


The Boltzman operator of eq. 
for the muon component is then 


B 


M* 




The neutron operator is 



(4.2) 

(14.1) 


Bn = n 


a 

dr + 


(14.2) 


The remainder are all simplified by the 
omission of charged particle stopping 


calculations to follow tnat tne power law 
model should be right on the average, 
than it should be right at some particular 
secondary energy or angle. 

The reaction cross sections used in 
the calculation to follow are assumed to 
be constant and geometric, i.e. a = nr§ l/A, 
ro=1.28 A3 f, where L is Avogadro's 
number. This assumption is valid for low 
energy nucleon -nucleus collisions for 
atomic weights from less than 12 to 
greater than 64 (ref. 12) , and it is used 
in what follows for all energies and 
hadrons. 


B = 0 + 

P dr 


(14.3) 


Bn = n L + CT + _£n_ (14>4) 

This omission is not important for 
secondaries above about 1 GeV (ref. 14) . 
Compensation for this can roughly be made 
with the use of the Heaviside function of 
eq. (7) , as will be shown. 


Separating the primary nucleons from 
the secondaries produced in the atmosphere 


^q = ^iq + cpsq (15) 
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Table III. - Shower particle multiplicities. 


E 

(GeV) 


Power Law 17 
model 

Statistical® 

model 

Extrapolation® 

model 

TRB^ 

model 

CKP 5 

model 

12. 5 

3.6 

5.0 

3.0 

5.5 

3.8 

4. 9 

20 

5.3 

5.9 

3.8 

- 

4.3 

5.5 

30 

7.3 

6.8 

4.8 

- 

5.0 

6.1 

200 

12 

12 

- 

6.9 

10 

9. 8 

300 

13 

13 

8.5 

" 

12 

11 

1000 

17 

17 

- 

" 

20 

15 


a nuclear emulsion. 

8 hydrogen, total charged particle yield. 

C aluminum, total charged particle yield, kaons neglected. 
& air, kaons neglected. 

S hydrogen, total charged pion yield. 


Where 

cp iq is the flux of primary nucleons, 

1 and 

cp sg is the flux of secondary nucleons. 
This leads to 3 differential equations 


[° 


— + — ti- - k u 

dr P a r dE ^ 


C n ± 

P~Tr CPsn± 

TT* 


1 "Vt* 

>sn = 

nVr± 


('• 


™\1 


2 17 

E 




n ) 


(16.1) 


n + a Ico = 0 

u or j j.q 


(q = p, n) 


(16.2) 



Equation (16.1) is written in integral 
form and reduced to quadratures (ref. 15) . 

A 51 point set was found necessary for the 
integral over space and 7 for the integral 
over angle. 

The solution to eq. (16.2) is 

cp iq = exp (- r/an) (17.1) 

Passow (ref. 16) , and Alsmiller (ref. 17) 
have shown that the solution to an integro- 


differential equation of the form of equa- 
tion (16.3) with constant cross section 
for incident nucleon flux of energy Eo and 
zenith direction Cl is 

■a 3 ( ^ \ 

’ TTT^ST) < l ' W ' 

' Pq rJ 

r r n* < 17 ’ 2 > 

• U(E ° - V L OB (Eo , E) -1 • 

• I 1 [2/ (r/0) B(Eo , E) ] 

B(Eq,E) = ff £ (1 - ^)Kt • 

t=n,p 

• {in E 0 - ln[EU(E - n t ) + Ti t U(Ti t - E) ]}, 
(v = p,n: q = p,n,n) 

Equation (17.2) differs slightly 
from the form obtained by Passow (ref. 16) 
and Alsmiller (ref. 17) by the inclusion 
of the decay term in the solution. The 
reason for this lies in the fact that 
C n = C p = 0, that G n j does not appear under 
the integral of eq. (16.3), and that pion 
production and absorption is purely local 
so that the energy independence of the 
cross section required by the solution can 
be relaxed. 

The remaining parameter, T|q, is a 
lower energy limit beneath which secondary 
particles of type q are suppressed. Thus, 
T) n = Eo , and r\ n = 0. To compensate for 
the neglect of proton stopping, rip is set 
equal to 500 MeV. 

The neglect of charged particle 
stopping, the straight ahead approximation, 
and the constant geometric cross section 
make eq. (17) increasingly shaky as 
secondary particle energies go below 1 GeV 
and fails completely by 100 MeV. At high 
energies, eq. (17) should become and 
remain quite accurate as long as l and a 
can be treated as constants. It may seem 
that the neglect of hadron production by 
incident pions should fail at high energy 
since (eq* 3) can become very long. 

However, as argued by Adair (ref. 18) , any 
pion Which interacts with a nucleus can be 
treated as an absorption. This arises 
because of the combination of the rela- 
tively low pion inelasticity with the 
steepness of the nucleonic energy spectrum. 
Only rarely will a pion be emitted from a 
pion -nucleus collision with an energy near 
to that of the incident nucleon (see Fig. 

1 for the predictions of CKP and TRB at 
high secondary energies) . The steepness 
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of the energy distribution then causes the 
number of pions resulting from pion- 
nucleus collisions to be small conpared 
with those produced directly in nucleon- 
nucleus collisions. 

The restrictions on the form of Gqj 
imposed by the power law model would seem 
a. Priori to be the most serious limitation 
on the use of eq. (17) . It is then worth 
pointing out that eq, (17) , suitably 
modified, has been applied to accelerator 
beam measurements in iron for proton 
energies from 1 to 18 GeV (refs. 5, 19), 
and an 8 GeV pion beam on tin (ref. 20) 
with excellent results. 

Once it is established that the param- 
eters underlying the power law model are 
correct and that the cross sections are 
correct, more general methods can be used 
and many of the approximations made here 
to stay within the bounds of Passow's 
mathematical framework (refs. 16, 17) can 
be abandoned. The relative ease and 
transparency of Passow's approximation 
make it of value in itself however. 

Photons and Electrons 

Electromagnetic shower propagation is 
not so problematic as nucleon transport. 
Essentially exact Monte-Carlo treatments 
exist and have been tested against experi- 
mental data (ref. 21, 22) . These calcula- 
tions are difficult to carry out over the 
range of depths and energies required, and 
as the goal at this time is the establish- 
ment of sufficient conditions to determine 
the atmospheric flux, the propagation of 
the electromagnetic cascade is treated very 
primitively. 

Since the mean life Q of the neutral 
pion is 0.91 x 10 ” 16 s, SlL_ is huge 
compared with a (c ^ 0 = 3.3 x 10 10 GeV) it 
decays immediately into 2 photons. The 
muon decay probability is very much less 
(C^ = 1.1 GeV) and is of significance only 
below about 10 GeV. Energy deposition is 
calculated from the assumption that the 
total energy of the neutral pion produced 
per gram of air (from eq. 17.2) is absorbed 
at the point of production, and 1/3 the 
total energy of the decaying muon (from 
eq. 16.1) is absorbed at the point of 
decay.* 


♦However, at this time, an attempt is being 
made to apply the electromagnetic shower 
code CASCADE (ref. 22) to this problem to 
improve the treatment of this important com- 
ponent of atmospheric cosmic rays. 


This assumption will deteriorate in- 
versely as geomagnetic latitude. As the 
geomagnetic cutoff rises toward 17 GV, the 
neutral pion production spectrum will 
become harder, and as the radiation length 
of air is of the same order of magnitude as 
the nucleonic collision mean free path, 
neglect of transport will become increas- 
ingly serious. However at higher latitudes, 
the error will be seen to be tolerable. 

SOLAR AND TELLURIAN MODIFICATION OF THE 
GALACTIC COSMIC-RAY SPECTRUM 

Solar Activity and the Interplanetary 
Medium 

It is well known that variations of 
solar activity, through the agency of the 
solar wind modulates the cosmic-ray 
spectrum found at the earth's orbit. This 
modulation is a consequence of cosmic-ray 
transport through the interplanetary medium 
and it is formally the same as that which 
would be produced by a heliocentric electric 

field having a magnitude at the earth's 
orbit of about 100 MV at solar minimum and 
about 1000 MV at solar maximum (ref. 23, 

24). 

The electric field model is a useful 
Computational tool for representing solar 
effects as the potential is the only 
adjustable parameter. This representation 
is for convenience only. It is not 
asserted here that a heliocentric poten- 
tial of this size exists and is responsible 
for solar modulation. 

The electric field model of the modu- 
lated cosmic-ray flux is (ref. 25) 



T = E + ZU 

P (x) = i Vx 8 + 2A nip c 8 x, 

W(x) = x + A mp.c 8 

vdiere 

no is the unmodulated galactic 

spectrum of atomic weight A, 
and atomic number Z, per 
steradian per era 8 per s per MeV, 
having an energy of E MeV, and 
U is the solar potential in MV. 

For the calculations to follow, the 
unmodulated spectrum is taken from Freier 
and Waddington (refs. 24, 26) below 10 GeV 
per nucleon. Above that energy the 
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spectrum is taken from Peters (ref. 27) 
with which it has a smooth overlap. 

Calculations of atmospheric ioniza- 
tion will be made and compared with some 
of the measurements performed by Neher 
(ref. 28) . These measurements have been 
analyzed to yield the incident proton 
cosmic-ray spectrum. This is compared to 
the predictions of the electric field 
model for U = 200 MeV in Fig. 3. The 

measurements were performed during the 
course of several months in 1965, and are 
near a solar minimum, but a minimum not as 
deep as that of 1954 (ref. 28) , so the 
value of U is reasonable. The agreement is 
quite good, and this spectrum will be used 
for the ionization calculations to follow. 
Some neutron calculations will require 
different values of U. These will be taken, 
where possible, from Freier and Waddington 
(ref. 26) . 



Fig. 3. The East-West variation of geomagnetic 
cutoff rigidity as a function of geomagnetic 
latitude cong>ared with the effective cutoff 
rigidity for an isotropic detector. 


The Earth's Magnetic Field 

The magnetic field of the earth 
deflects incoming cosmic-rays depending on 
their rigidity and angle of incidence, so 
that for each angle of incidence there is 
a critical rigidity below which the incom- 
ing particle cannot interact with the 
earth's atmosphere. In the calculations 
which follow, a single cutoff rigidity will 
be applied. The primary spectrum will be 
assumed unchanged in angle and energy 
above the cutoff, and vanish below it. 
Richtmeyer et al. (ref. 29) have calculated 
the effective cutoff rigidity seen by an 
isotropic detector exposed to the primary 
spectrum at the top of the atmosphere and 
this cutoff will be used here. 


In Fig. 4 the cutoffs for the eastern 
and western horizons and the zenith, from 
Lemaitre and Vallarta (ref. 30) are shown 
along with the values obtained by 
Richtmeyer et al. (ref. 29) as a function 
of geomagnetic latitude. It is evident 
that the assumption of isotropy of the 
radiation near cutoff is not justified at 
latitudes below 40° to 45° . 

This, in combination with the assump- 
tions made with respect to electromagnetic 
shower transport, will probably cause the 
calculations based on an isotropically 
incident spectrum to fail at low latitudes. 
Consequently, at this stage in the develop- 
ment of the calculations the experimental 
compar isons will be limited to higher 
latitudes. 



Fig. 4. comparison of electric field model , 
calculations of the incident coBmic-ray proton 
spectrum with Usher's measurements. 


COMPARISON WITH EXPERIMENT 

Sea-Level Particle Spectra 

The first step in producing integral 
quantities such as ionization or neutron 
density is the calculation of the differ- 
ential energy angle particle distributions. 

As a test of eq. (17) , the .vertical 
component of the cosmic-ray spectra has 
been calculated for a geomagnetic 
latitude of 57°. Equation (17) with 0=1 
and r = 1033 g/cm a , was integrated over 
the source spectrum. Seventy percent of 
the source spectrum was assumed to be 
composed of free protons, and 30 percent 
of bound neutrons and protons all having 
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the energy distribution given by eq. (18) , 
with U - 200 MV. For geomagnetic purposes 
bound nuclei were treated as bound, but 
for the purpose of atmospheric transport 
treated as free i.e., an a particle is 
assumed to behave exactly like 2 free 
neutrons and 2 free protons. In Fig. 5, 
the calculated vertical component of the 
cosmic-ray nucleon spectrum of one charge 
state (neutrons or protons) is compared 
with the experimental sea level proton 
spectrum of Brooke and Wolfendale (ref. 31) 
and the sea level neutron spectrum of 
Ashton and Coats (ref. 32) . 

At this atmospheric depth, eq. (17) 
predicts nearly equal numbers of neutrons 
and protons, and so both experimental and 
theoretical data were combined. Agreement 
is very good over 4 decades of energy and 
10 of intensity. 

The calculated vertical component of 
the sea level pion spectrum is compared 
in Fig. 6 with the measurements of Brooke 
et al. (ref. 33) . Agreement is good over 
most of the range of comparison. 

The sea level muon spectrum for a 
zenith angle of 0° is shown in Fig. 7 
compared with measurements by Owen and 
Wilson (ref. 34), Holmes et al. (ref. 35), 
Gardener et al. (ref. 36) and Hyman and 
Wolfendale (ref . 37) . Agreement with 
experiment is quite satisfactory over the 
range from about 1 to about 1000 GeV. All 
measurements and calculations were for 57° 
geomagnetic latitude. 



)©■« >0® to' 10* 10* 10® 

K««Vt 

Pig. 5 . The vertical conponent o£ the cosmic-rey 
nucleon flux of one charge state at sea level as 
measured by Brooke and Wolfendale (ref. 31 ) , by ' 
Aahton and Coata (ref. 32 ), and as calculated. 



E (GeV) 


Fig. 6. The vertical conf)onent of the cosmic-ray 
pion flux as measured by Brooke et al. (ref. 33), 
and as calculated. 



Fig. 7. The vertical component of the cosmic-ray 
muon flux as measured by Owen and Wilson (ref. 

34) , Holmes et al. (ref. 35), Gardener et al. 
(ref. 3 6) , and Hayman and Wolfendale (ref. 37) , 
and as calculated. 
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In Pig. 8, the sea-level muon spectrum 
is calculated for a zenith angle of 75° and 
compared with the measurements of Stefanski 
et al. (ref. 38) . The lowest experimental 
point is above the maximum cutoff for this 
latitude, 52° (see Fig. 4) . 

As has been observed earlier, the 
assumptions that lead to eq. (17) become 
increasingly shaky below about 1 GeV. 

This affects all the charged particle 
distribution, including the muon energy 
distribution which of course depends on the 
pion distribution, and can be seen clearly 
in Fig. 5, 6, and 7, vhere the calculated 
fluxes rise above the measured fluxes in 
every case. 



Cosmic -Ray Ionization in the Atmosphere 

Ionization from protons, pions and 
iraions is calculated by multiplying the 
energy distributions by the appropriate 
stopping powers as described earlier (ref. 
1) . The lower energy limit for the proton, 
charged pion and neutral pion energy dis- 
tributions is 100 MeV, below vftiich the 
theory fails. The muons are allowed to 
slow down to 10 MeV, below which very 
little is contributed to the ionization. 

It was found necessary empirically tc 
use an upper limit of 10 4 GeV to include 
all significant contributors to the 
ionization. 


In Fig. 9, the calculated ionization 
(in units of I, the number of ion pairs per 
cm 3 of NTP air) at a geomagnetic latitude 
of 55° is compared with the measurements of 
Neher (ref. 28) , and later data as reported 
by George (ref. 39) down to 600 g/cm 3 . To 
complete the curve, the results of Lowder 
and Beck (ref. 40) from 600 g/cm a to sea 
level measured at 51° geomagnetic latitude 

at about the same time are included. 



Over-all agreement is seen to be 
within 20% with the exception of the 
region near 600 g/cm 3 , where the disagree- 
ment is nearer 40%. The comparison is 
absolute it must be emphasized. The 
composition of the total ionization is 
shown in Fig. 10. Because the secondary 
fluxes interact differently with the 
atmosphere each component has a noticeably 
different profile. The kink in the 
electron curve about 850 g/cm 3 is a con- 
sequence of the transition from shower 
production originating in neutral pion 
decay at low depths to shower production 
resulting from muon decay at larger depths. 
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A similar comparison is shown in Fig. 
11, where the calculation is carried out 
at 44° , and the data again are taken from 
Neher (ref. 28) . To complete the curve, 
the data of George (ref. 39) are included 
from 188 g/cm 8 to sea level. The disa- 
greement is typically 20% with higher 
values at near 800 g/cm 8 and 50 g/cm s . The 
latter is probably a result of the depar- 
ture of the incident primary cosmic-ray 
flux from isotropy, and the hardening of 
the photon production spectrum that 
results from the higher average cutoff 
(see Fig. 4). George's (ref. 39) measure- 
ments were carried out during 1968 near a 
solar maximum, and hence the additional 
modulation, if removed, would make the 
disagreement worse. 



Pig. 11. The cosraic-ray ionization profile at 
44 as measured by Neher (ref. 28) and by George 
(ref. 39 ) and as calculated. 


Cosmic-Ray Neutrons in the Atmosphere 

The lower energy limit of the cascade 
calculations described above is 100 MeV, 
and the agreement with measurement 
indicates that this cutoff which is forced 
on the calculations by the limitations of 
the analytical theory, is adequately high. 
This is a consequence of charged particle 
stopping which limits the number of 
charged particles at low energies. Neu- 
trons are uncharged however, and extend 
all the way down to thermal energies. In 
order to account for neutron fluxes below 
100 MeV, the cosmic-ray neutron spectrum 
reported by Hess et al. (ref. 41) at sea 
level and 44° geomagnetic latitude has 
been patched onto the calculated differen- 
tial spectrum at 100 MeV. This approach 
is rather rude, and fails at small depths 
as it cannot account for the diffusion 
hardening which takes place near a vacuum 
boundary. 


This is clearly seen in Fig. 12, where 
the neutron flux measurements of Boella et 
al. (refs. 42, 43) , and Yamashita et al. 
(ref. 44) are compared with calculations. 

At depths greater than 200 g/cm s the 
agreement is really rather good. The lack 
of isotropy of the incident flux near cut- 
off and the absence of diffusion hardening 
lead to an overestimate at small depths. 

The measurements of Boella et al. (ref. 42) 
and Yamashita et al. (ref. 44) are both 
ground level measurements rather than free 
air measurements, but yet are seen to fall 
on the curve. The effect of the air ground 
interface on the calculations has not been 
evaluated. 



Pig. 12. Ww co*mic-ray neutron flux in the 
atmosphere as measured by Boella et al. (ref. 

42) , Boella et al. (ref. 43) , and Yamashita 
et al. (ref. 44), and an calculated. 

In Fig. 13, the neutron density meas- 
urements of Yuan (ref. 45) and Gold (ref. 

46) have been compared with calculations 
on the same basis. Again, the agreement 
between calculations and measurements is 
good except at small depths, with the 
exception of the ground level value of Gold 
(ref. 46) . Gold (ref. 46) recognizing that 
his values were quite high compared with 
balloon measurements in free air (ref. 47) 
attributed this to the interface with the 
ground, to which the neutron density 
appears much more sensitive than the neu- 
tron flux. 

Yuan 1 s data appear to have a markedly 
different slope from the measurements. 

This can also be seen in Boella et al. 

(ref. 43) on Fig. 12. This effect may be 
a consequence of operating near the 
threshold of instrument sensitivity since 
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the calculation is in good agreement with 
the deeper measurements, all the way down 
to sea level both in Fig. 5 and 12. Such 
an effect will lead to rather long 
reported relaxation lengths. 



Fig. 13. The cosmic-ray neutron density in the 
atmosphere as measured by Yuan (ref. 45) , and 
Gold (ref. 46) , and na calculated. 


Figure 14 shows the data of Miles 
(ref. 48) . In this case a general under- 
estimate appears below 200 g/cm 2 . The 
identical calculations can be applied to 
the data of Haymes (ref. 49) and Soberman 
(ref. 50) which have been converted to 
neutron density and adjusted by Miles 
(ref. 48) for differences in latitude and 
time so that comparison with his own re- 
sults might be made. Above 100 g/cm a , the 
data of Soberman (ref. 50) in Fig. 15 are 
seen to be in excellent agreement with 
calculation. At smaller depths, the 
neglect of leakage, which was not taken 
into account in the low energy model, 
causes the calculation to be too high. 
Although much more scatter appears in the 
measurements of Haymes (ref. 49) , agree- 
ment is clearly reasonably good except 
again at small depths. 

DISCUSSION 

Analytical calculations of the 
secondary energy distribution of cosmic- 
rays in the atmosphere, cosmic -ray ioniza- 
tion, and neutron flux and density have 
been performed. The model of nucleon- 
nucleus collisions on which the calcula- 
tions depend is based on the constant 
partial inelasticities of Table I, and 
constant geometric reaction cross sections. 
The transport calculations assume a purely 
nucleonic cascade, and consider (eq. 5) 
only protons, neutrons, pions, electrons 
and photons. 



Fig. 14. The cosmic-ray neutron density in the 
atmosphere as measured by Miles (ref. 40), and 
as calculated. 



Fig. 15. The cosmic-ray neutron density in the 
atmosphere as measured by Soberman (ref. 50) , 
and adapted by Miles (ref. 48) , and as 
calcul ated. 

In particular; the muon ionization 
and the high zenith angle energy distri- 
butions calculated on the basis of tt* -* [l 
decay, are in agreement with measurement 
(Fig. 8) and support the conclusion of 
Stefanski et al. (ref. 38) that there is 
no other important contribution to the 
muon flux below 300 GeV such as would 
account for the Utah muon measurements 
(ref. 51) . 
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Ranft and Borak (ref. 9) observe that 
the measurements of negative pions at 70 
GeV of Bushnin et al. (ref. 52) are in 
strong disagreement with the TRB model of 
hadron nucleus collisions and with the 
statistical model of strong interactions 
(ref. 3) . The extrapolation model disagrees 
by an order of magnitude with the same data 
(ref. 8). If the measurements are correct, 
they cannot imply a change in the appropri- 
ate partial inelasticities. As has been 
shown the power law model used in the 
calculations of the fluxes in figs. 5-8 is 
in agreement with these models. Thus a 
large error in the assigned values of Vj or 
Kj would lead to quite large errors in the 
differential nucleon, pion and muon fluxes 
starting somewhere between 30 and 70 GeV. 
However, the measurements are only at small 
angles (^ 15 milliradians) and high sec- 
ondary momenta (> 45 GeV/c) , and it is 
possible that average quantities such as 
multiplicity and inelasticity are not 
much affected by what happens in this 
region. 

It has been suggested that nucleon- 
nucleon and meson-nucleon cross sections 
may vanish at infinite energy (ref. 53) . 
These calculations, performed with constant 
cross sections, indicate that hadron-nucle- 
us cross sections are essentially constant 
and geometric out to 10 4 GeV. For instance, 
a 10% perturbation of the nucleon-nucleus 
cross section will cause a change in the 
nucleonic flux in the region of 10 3 - 10 4 
GeV of 150%. This would put theory and 
experiment out of agreement in Fig. 5. 

But, such a discrepancy could be accounted 
for in terms of errors in the primary 
spectrum used, or errors in the sea level 
measurements. Much larger changes in the 
cross section however, would lead to very 
big changes in the sea-level flux which 
would not be reconcilable with the data. 

In, these calculations, all hadron- 
nucleus cross sections are equal and 
geometric. It is known that pion-nucleon 
cross sections are 2/3 of the nucleon- 
nucleon cross sections (ref. 53). As 
Adair (ref. 18) points out, pion-nucleus 
and nucleon-nucleus cross sections will 
differ by less than this as a consequence 
of the intra-nuclear cascade. He calcu- 
lates the ratio to be 0.77. Ranft and 
Borak (ref. 9) obtain 0.83 from published 
experimental data. The effect on pion and 
muon spectra of a 20% error in the reac- 
tion cross section would be much less than 
in the nucleon case, as pions are locally 
produced, and locally absorbed. 


Apparently one-dimensional nucleonic 
cascades depend only weakly on the angular 
and energy behavior of the secondary 
production spectra, but strongly on the 
partial inelasticities, multiplicities and 
cross sections. The calculation described 
here which is based on a straight ahead 
power law approximation to the doubly 
differential production spectrum yields 
good agreement with the sea level differ- 
ential fluxes of nucleons, pions and muons, 
with the cosmic-ray ionization at all 
depths in the atmosphere, and with the 
neutron fluxes and densities in the 
atmosphere. 

Lastly, it may be observed that in 
Figs. 12, 13, and 14, the neutron calcula- 
tions appear to relax more rapidly than 
some sets of measurements. Yet in the 
Fig. 5, and the sea level values of Fig. 

12, this, would appear to represent 
experimental error, possibly a consequence 
of operating near the threshold of experi- 
mental sensitivity. Too rapid a relaxa- 
tion rate of the nucleon spectrum would 
also result in the calculated pion 
spectrum of Fig. 6 to be too low, which it 
is not, as the pions are produced locally 
from nucleon-nucleus collisions. 

The neutron attenuation is not 
actually exponential although its 
departure from exponentiation is not 
large (ref. 54) . In Table IV, relaxation 
lengths obtained from the neutron density 
calculations at X = 41° and with U = 400 MV 
are compared with Miles (ref. 47) . Agree- 
ment is poor, but the calculated neutron 
densities agree well with Haymes (ref. 49) 
and Soberman (ref. 50) when they are 
reduced to the same conditions and are 
never less than half of Miles (ref. 49) . 

In addition, Simionati deFritz and Cicchini 
(ref. 54) have measured cosmic-ray neutron 
attenuation lengths in air at 25° (ref. 55) 
as a function of atmospheric depth which 
are included for comparison. In this case 
agreement is quite close, and this suggests 
that exponential relaxation lengths are not 
well determined by experimental attenuation 
data, as relaxation rates are not constant 
with height and the true variation is 
marred by poorer quality at greater depths 
and at lower intensities. 
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The success of the preceding calcula- 
tions rests primarily on the phenomenolog- 
ical model of hadron- nucleus collisions 
presented here. The description of the 
exclusively nucleonic cascade combined with 
the neglect of kaons is a simplification 
with which more sophisticated and expensive 
calculations may dispense. Other, more 
obvious, approximations are being discarded 
as the study progresses. 


Real improvement on the forms of the 
doubly differential production spectra 
must probably await progress in studies 
being carried out at Oak Ridge and else- 
where. The adequacy of the power law 
nuclear model for calculations of this type 
is clear from the results shown here, and 
the average properties, particularly the 
partial inelasticities and cross sections, 
as applied to air, appear to be well 
established. 


Table IV. - Comparison of calculated neutron 
attenuation lengths for X = 41° and U = 400 
MV with measurements. 


Atmospheric depth 
(g/cm 3 ) 

Calculated 

(q/cm 8 ) 

Measured 

(g/cm 3 ) 

300 

149 

155<7 

500 

129 

125 a , 165±20£ 

1033 

113 

115 a 


a X = 25° (ref. 54) . 
3 X = 41° (ref. 47). 



Fig. 16. The cosmic-ray neutron density in the 
atmosphere as measured by Haymes (ref. 49) and 
adapted by Miles (ref. 48) , and as calculated. 
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